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Abstract. We present arguments in favour of the inequalities var(X^ | X S B v (p)) < 2A n E[X^ | X S B v (p)], where 
X ~ Afv(0, A) is a normal vector in V > 1 dimensions, with zero mean and covariance matrix A = diag(A), and B v (p) is a 
centered v— dimensional Euclidean ball of square radius p. Such relations lie at the heart of an iterative algorithm, proposed 
in ref. [1] to perform a reconstruction of A from the covariance matrix of X conditioned to B v (p). In the regime of strong 
truncation, i.e. for p < A n , the above inequality is easily proved, whereas it becomes harder for p> A„. Here, we expand both 
sides in a function series controlled by powers of X n /p and show that the coefficient functions of the series fulfill the inequality 
order by order if p is sufficiently large. The intermediate region remains at present an open challenge. 



1 Introduction 

It is intuitively clear that independent random variables develop correlations once constrained within compact 
multivariate domains. Whenever the mathematical framework rules out closed-form results, a possible 
approach to studying such correlations is to focus on inequalities among expected values. As a case in 
point, in this paper we consider a random vector X ~ M v {0, A) in v > 1 dimensions, with A = diag(A) and 
A = {Afc}^ =1 , whose probability density is cut off sharply outside a Euclidean ball 

B v (p) = {xeR v : x T x < p}. (1.1) 

Owing to the symmetry mismatch between Af v (0, A) and B v (p), the conditional moments of X admit no 
exact representation in terms of elementary functions. Our aim is to show that the effect of the spherical 
truncation on the variance of the square components of X is quantified by the inequalities 

A n (p; A) = ~ {var (X 2 n \ X G B v (p)) - 2X n E[X 2 | X G B v (p)]} < , n = 1, . . .v. (1.2) 



The interest we have in eq. (1.2) originates from ref. [T], where we have proposed a fixed-point algorithm for 



the reconstruction of A, in case the only available information amounts to the covariance matrix &b of X 



conditioned to B v (p). In particular, in that paper we showed that eq. (1.2 1 and the correlation inequalities 
cov(Xl X 2 m | X e B v (p)) < , n £ m , (1.3) 
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Fig. 1 - Contour plots of A„ at v = 2. 




are necessary and sufficient for the convergence of the algorithm. Eq. (1.3) expresses a property of negative 
association among the square components of X (see ref. [2]). A proof of it goes beyond the scope of the 
present paper. 

If we denote Gaussian integrals over Euclidean balls by 



otkim...{P\ A) 



,y>2 ,-y, 2 

d v x — — 

Afc A^ A m 



, A,) 



e -y 2 /(iv) 
(2tt?7) 1 /2 



(1.4) 



and define <9„ = d/(d\ n ), we see that A„ = 2(A 2 l / / o 2 )[A„d n (a„/a)]. Thus, the inequality A„ < holds true 
iff a n /a is monotonic decreasing in A„ with A(„) = {Ai}j^ n kept fixed. Since such monotonic behavior is 
held by both a n and a separately, eq. ( |1.2[ ) simply means that a„ is more rapidly decreasing than a. An 
illustrative example is shown in Fig. [l] where contour plots of A n at v — 2 are reported^] 



Apart from the covariance reconstruction algorithm, a different motivation to care about eqs. (1.2) and 



( 1.3 ) has to do with non-linear optimization issues. Thanks to Prekopa's Theorem [3J, a(p; A) is easily shown 
to be logarithmic concave in p. In sect. 2 we discuss how log-concavity relates to the correlation inequalities. 



Now, despite the seemingly candid aspect of eq. ( |1.2 
which works across the whole parameter space (p; A) G 



The outcome is that, independently of Prekopa's Theorem, eqs. (1.2) and (1.3) are alone sufficient to induce 
log-concavity, while they cannot be deduced from it. 

it is difficult to find a unique rigorous proof of it, 
x . In this paper we propose three different 
partial arguments. The first two are discussed in sect. 3. They are both straightforward and apply in the 
regime of strong truncation, i.e. for < p < A n resp. < p < 2A„, independently of A( n ). In particular, 
the first one is based on Holder's inequality, while the second one makes use of the integral representation of 
A n . The relative ease of proving eq. (1.2) at strong truncation is certainly due to the large negative values 



A n assumes in this regime and its weak dependence upon A(„j, which in a sense makes the problem nearly 
1 -dimensional. 



The third argument, presented in sects. 4 and 5, applies instead in the regime of weak truncation, i.e. for 
p >• A„, where proving eq. ( 1.2 ) is definitely harder. As p — » o o, we have indeed var(A 2 | X e B v (p)) — > 2\\ 
and E[A, 2 j | X g B v (p)] — > A„, thus A„ — > 0. Hence, if eq. (1.2) is correct, it must follow from a cancellation 
of two positive terms resulting in an increasingly small negative balance. Motivated by the observation 
that the volume constraint weakens as p — > oo (and consequently akim— becomes well approximated by a 
product of 1-dimensional Gaussian integrals), we expand A„ in a non-elementary-function series around 
the factorization point. Each term of the expansion factorizes into a 1-dimensional integral along the n th 
direction plus a residual (v— 1) -dimensional integral in the orthogonal subspace. We prove that such factors 
get opposite signs as p — » oo, thus resulting in negative contributions. 

1 Numerical techniques for the computation of akim... are discussed in ref. pQ. 
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Wc finally draw our conclusions in sect. 6. 



2 Basic properties of OLMm... 

It is worthwhile starting our discussion by reviewing some trivial properties of the Gaussian integrals, which 
are used in the sequel. As an alternative notation for akim... we sometimes adopt the symbol a.\;ki...v.k v i 
where kj counts the multiplicity of the index j = 1, . . . , v. Whenever a directional index has zero multiplicity, 
we simply drop it. For instance, we write ay : fc. in place of the more pedantic ai:o...j:kj...v.o- When needed, 
we declare the integral dimension of ctkim... explicitly by writing the latter as a ktm . With this in mind, we 
proceed to a first set of statements: 

Proposition 2.1. Gaussian integrals fulfill the following properties: 

Pi) aki m ...{p]X) is increasing in p; 
P2) ctktm...{Pi X) * s separately decreasing in Ai, . . . , X v ; 
P3) a kim... {p'i A) fulfills the scaling equation 



pd p + ^ Kd r 



r=l 



a Hm ...(p;X) = 0; (2.1) 



p^) one-index integrals ctk-.n follow the hierarchy 

a k .. n < (2n - l)a fc! ( w _i) < (2n - l)(2n - 3)a fe:( „_ 2) < . . . < (2n - 1)!! a < [2n - 1)!! ; (2.2) 

P5) differentiating aki m ... (p; A) with respect to p yields 

1 1 " 

pd p a kl ... kn {p; A) = -{v + 2n)a kl ...k n (p; A) - - ^ ^...k^ip; A) , n = 0, 1,2, ... (2.3) 

fe=i 

Pq) a{p; A) is logarithmic concave in p, i.e. it fulfills 

a(s Pl + (l- s)p 2 ;X) > [a(p 1 ;X)] s [a(p 2 ;X)] 1 ~ s , 0<s<l, Pl ,p 2 eR+. (2.4) 

Proof. Property pi) follows from the positiveness of the integrand of akim... and the observation that if 
Pi < Pi then B v (pi) C B v (p 2 ). H Property p 2 ) follows from the observation that akim... depends on p and 
A only via adimensional ratios, as seen by rescaling the integration variable x — > xj ^fp in eq. (1.4), i.e. 



( u P V/2 f rv P^tpxjpXm j Pf^l 

otum...{p, a) = (27r) „ /2|A|1/2 y Bt w d x— — — . . . ex P j - - x; ^ 1 

= f 1 ; {7. -.7}) • ( 2 - 5 ) 

When a single variance is downscaled, e.g. A — > A' = {Ai, . . . , aX r , . . . , A^} with < a < 1, the change in 
aktm... is entirely transferred to the integration region, i.e. 

aMw ~ (p; A ) = (W d x V V XT ' ' ' exp r 2 £ ) - (2 - 6) 
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with 



£ v (l; a) = {x £ R v : x\ -\ h ax 2 r -\ h^<l} 



(2.7) 



Since B v (l) C £ v (l;a), it follows otklm...{p'i A') > ay ra ...(p; A). □ Property ^3) follows from the application 
of the chain rule of differentiation to eq. (2.5). The meaning of the scaling equation is that akim... keeps 



invariant under a change of the units adopted to measure both p and A. □ To get convinced about property 
Pi), we first notice that 



Afc^fe < ^l:ni...k:n^...v.' 



1 



[&l: ni ...k:(n h +l)...v:n v ~ {^ n k + l)ai:m ■■■k:r 



(2.8) 



as proved by evaluating the derivative on the l.h.s. under the integral sign. From property p 2 ), the r.h.s. of 
eq. (2.8) is recognized to be negative. The proof is completed by taking rij — for j 7^ k. Note that eq. (2.2) 



entails the inequalities 

E[Xf |XeS»] < (2n-l)!!A 



k 1 



n=l,2,. 



(2.9) 



with the quantity on the r.h.s. representing the value of the unconditioned (2n) th univariate moment of X^. 
We shall use the lowest order inequalities E[X% \ X £ B v (p)] < A& and E[X% \ X £ B v (p)] < 3A^ time and 
again in the sequel. Note also that the larger n, the slower otk-.n saturates to its infinite volume limit. Indeed, 
if we denote by d n = [(2n — 1)!! — otk-.n]/ Q-k-.n the fractional distance of otk-.n from its infinite volume limit, 



then eq. (2.2) is equivalent to the inequality chain 



d (p; A) < d x ( P ; A) < d 2 ( P] A) < . . . . (2.10) 
As we shall see, this property lies at the heart of most of the difficulties related to proving eq. (KL2j) . □ 



Property p§) follows from eqs. (2.1) and (2.8). □ Finally, in order to prove property pe), we recall 



Theorem 2.1 (Prckopa). Let Q(x) be a convex function defined on the entire v -dimensional 
space W" . Suppose that Q(x) > a, where a is some real number. Let ip(z) be a function 
defined on the infinite interval [a, 00). Suppose that iJj(z) is non-negative, non-increasing, 
differentiable, and —ip'(z) is logarithmic concave. Consider the function f{x) = tJj(Q(x)) 
(x S W ) and suppose that it is a probability density, i.e. 



d v x f{x) = 1 



(2.11) 



Denote by P{C} the integral of f{x) over the measurable subset C ofW. v . If A and B are 
any two convex sets in W , then the following inequality holds: 



(P{A}) S (PiB}) 1 - 8 < P{sA + (1 - s)B} , < s < 1 , 
where the linear combination on the l.h.s. denotes the Minkowski sum 
sA + (1 - s)B = {sx + (1 - s)y : x e A, y e B) . 



(2.12) 



(2.13) 



Obviously, theorem 2.1 applies if f(x) is a product of univariate Gaussian densities, as is the case with 



a(p;\). In addition, if x £ B v (pi) and y £ B v (p2), from the convexity of the square function x h-> x it 
follows that 



k=i 



sx k 



(1 - s)y k ] 2 < s J2x 2 k + (l-s)J2yl<spi + (l- s)p 2 



(2.14) 



k=l 



fc=l 
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i.e. sB v (pi) + (1 — s)B v (p2) Q B v {sp\ + (1 — s)/^)- Accordingly, we conclude that 



[a(p 1 ;X)] s Hp 2 ;X)] 1 - s < / d v x TT 5(x m ,X m ) 

/ s e„( Pl )+(i-6')e„(p 2 ) m=1 



< 



/ II 5 ( x m,X m ) = a(spx + (1 - s)p 2 ;X) ■ 

JB v {s P1 + {X-s) P2 ) m=1 



(2.15) 
□ 



Now, log-concavity is a local property of a(p; A), yet it brings global information about the conditional 



moments of X. To see this, we observe that since a(p; A) is twice differentiable with respect to p, eq. (2.4) 
is equivalent to 



adla 



(d P af<0. 



(2.16) 



We iterate eq. (2.3 1 to express the above derivatives in terms of conditional expectations. In first place, 



evaluating that equation at n = yields 

■ E[X*\XeB v (p)) 



2p 2p 



E 

k=l 



Xi 



(2.17) 



Property p\) then implies 
^E[Xi\XeB v (p)} 



k=l 



x k 



<v. 



(2.18) 



Though trivial, eq. ( 2.18[ ) calls for two remarks. The first one is that a sufficient (but not necessary) condition 
for it to hold true is E[Xl \X G B v (p)} < Xk Vfc, which has already been established. In second place, 
differentiating a in p an arbitrary number of times generates always symmetric expressions with respect to 
the directional indices, since p is not tied to any specific direction. In particular, this is the case with the 
second derivative, 



2 af v(v-2) E[X% | X £ B v (p)] 1 ^ E[A|X fc 2 | X € B v (p)] \ 

" a P 2 \ 4 2^ A, XjX k } 



(2.19) 



We see that all directional indices are again summed over. We shall come back in sect. 4 to the rational 



coefficients multiplying the expectation values on the r.h.s. of eqs. (2.17) and (2.19). For the time being, we 



finalize our argument by inserting these expressions into eq. (2.16). A little algebra yields 



1 k=l L j^k 



cov(XlXl\XeB v (p)) 



A, A 



< 0. 



(2.20) 



Eq. (2.20) describes the log-concavity of a in terms of conditional expectations. Did we not know about 



Prekopa's Theorem, we could regard it as a result of eqs. ( 1.2 ) and ( 1.3 ). Unfortunately, the converse does not 



hold: it is not possible to infer eqs. (1.2) and (1.3) from eq. (2.20), as contributions along different directions 



could compensate while keeping the l.h.s. negative. Nevertheless, if eqs. (1.2) and (1.3) were simultaneously 



violated for all indices, a(p; A) could not be logarithmic concave at all. Therefore, eq. (2.20) tells us that at 
least some of the correlation inequalities must hold. To conclude, a full proof of eq. (1.2) cannot follow from 
the property of log-concavity, so we need to look elsewhere. 
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3 Variance reduction in the regime of strong truncation 



For the sake of conciseness, throughout this section we denote conditional expectations by E[-] instead of 
E[ • \X € B v (p)]. Our starting point consists in regarding eq. (1.2) as an upper bound to E[X^]. This 
suggests to consider the wider inequality chain 



E[X% < E[X*} (2A„ + E[X*}) < X n (2A„ + E[X%) < 3Xl ■ 



(3.1) 



The leftmost bound is in fact a recast of eq. (1.2) 



If for a moment we give it for granted, the second and 
third ones follow as an immediate consequence of E[X^] < A„. Although our final target is just represented 
by eq. (1.2), it makes sense to first consider the two rightmost bounds: if they turn out to be violated, 
eq. (1.2) cannot be correct. The loosest one is once more the trivial inequality E[X„] < 3A^, which we have 



already established. By contrast, the inequality 
E[X*] < \ n (2\ n +E[X*]) 



(3.2) 



is less obvious. In sect. 3.1 we prove it. Our argument is based on straightforward algebraic manipulations 
of the Gaussian integrals over B v (p). We include it in the present note for a twofold reason: on the one hand 
it gives a feeling of the optimality of eq. (1.2), on the other it represents the only general result we have, 
valid across the whole parameter space. 



3.1 A loose yet general bound to E[X* | X e B v (p)} 



In order to prove eq. (3.2 ), we use a standard trick, consisting in a rescaling of A„ by an external parameter r, 

on 

(3.3) 



so as to obtain the moments of X n by differentiation of a in r. More precisely, we introduce the function 
1 



H(r) 



a\p 



Ai, A„l) = ^= [ d v x 8(x n ; X n /r) TT S(x m , A m ) , 
)) ^ JbAp) ™i 



whose dependence upon p and A we leave implicit. Differentiating H(t) under the integral sign yields 
,2 fc A£ d k H 
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E[Xf ] = (-1) 



a dr k 



k = 0,1,2,. 



(3.4) 



At the same time, derivatives of H{r) can be taken via the chain rule of differentiation, which allows us 
to express them as algebraic combinations of a and its derivatives in A„. For instance, with regard to the 
second and fourth moments, we find 



2r 3 / 2 



dU 
dr 

d 2 n 

dr 2 ~ 4r 5 / 



1 



(a + 2X n d n a) , 



^ (3a + 8X n d n a + 4X 2 n d 2 a) . 



(3.5) 



(3.6) 



Consider first the lowest order derivative. By inserting eq. (3.5) into eq. (3.4) evaluated at k = 1, we obtain 
E[A^] = A„[l + 2(X n /a)d n a\. Comparing this with E[X^] = X n (a n /a) yields 



a„ = a + 2X n d n a . 



(3.7) 



Eq. (3.7) coincides with eq. ( |2.8| ) evaluated at n\ = ... = n v = 0. Owing to property p 2 ) of sect. 2, we infer 
a n < a and thus we find again E[X„] = X n (a n /a) < A„. Consider then the fourth moment. If we insert 
eq. (3.6) into eq. (3.4) evaluated at k — 2, and then make use of eq. (3.7), we easily arrive at 



EK]=4A„E[^]-A^ + 4A^ 



4 d n a 



(3.8) 



In order to estimate d^a, we differentiate both sides of eq. (3.7 1 with respect to A„. We then invoke again 
property pi) of sect. 2, thus obtaining 



dl a 



1 3 

— — (d n a n - 3d n a) < -——d n a 



4A2 



(a n - a) 



3a 
iXj 



{E[X 2 n ] -A„} 



(3.9) 



This estimate is sufficient to prove eq. (3.2 1. 



3.2 First argument in favour of eq. (1.2) 



In the regime of strong truncation, eq. (1.2 1 can be inferred from Holder's inequality. We recall that if 



p, q > 1 are two numbers satisfying 1/p + 1/q = 1 and X , Y are stochastic variables on a given probability 
space, then E[\XY\] < (E[\X\p]) 1 /p(E[\Y\ c ']) 1 ^ . In our case, we have 

var(^) = E [(Xl E[X 2 n ]f] = E [X n ] E [A 2 ] 2 = E[X n - E[A 2 ] 2 ] 

= E [{X 2 n - E[X 2 n ]) {Xl + E[Xl])} < E [| {Xl - E[Xl]) (X 2 n + E[X*\) |] 



< 



{E [\Xl - E[Xl]\ p } } 1/P {E [{Xl+E[Xl]f] } 



1/9 



(3.10) 



The latter inequality holds true for any finite choice of p, q, provided their reciprocals sum to one. Accordingly, 



it holds as well in the joint limit q 
var(X 2 ) <2^E[A 2 ], 



with 



1 , p = q/(q — 1) — > oo, where it amounts to 



h = lim (E \\X 2 n - E[X^]| P 1 } 1/P = ess sup C\X 2 n ~ E\Xl] 



(3.11) 



(3.12) 



Recall that the essential supremum of a real-valued function / is defined by ess sup / = inf{a e K : [i{{x : 
f{x) > a}) = 0}. In particular, the measure /i which is understood in eq. (3.12| is the marginal probability 
measure of X n , i.e. 



, , > _ a ( " 1) {p-x 2 n ;X {n) )5{x n ,X n ) 
afi(x n ) — , . ax r , 

aS v >{p; X) 



(3.13) 



Owing to the modulating factor oS v ~ x \p — x^;A( n )), [i is neither Gaussian nor log-concave (in x n ). In 
sect. 4 we shall say more about eq. (3.13) and the factorization of its numerator into functions of resp. A„ 
and A(„). For the time being, we observe that [i has support in the interval (— ^/~p, +V7')- Depending on how 
E[X r 2 J compares with p/2, h might assume one of the values h\ = p — E[X^\ or hi = E[X 2 ], as qualitatively 
represented in Fig. [2] As far as we are concerned, we do not need to establish which among Fig. [2^, and 
Fig. [2Jd provides the correct qualitative behavior for 2 t — E[A r 2 J | : numerical computations suggest 

that Fig. [2}d is not realized for any choice of n, p and A, yet this information is irrelevant for what follows. 
More precisely, we distinguish three cases: 



i) P < A„ {strong truncation): in this case h < A„. Indeed, since E[A 2 ] < A„, both hi and hi lie below 
A n . In this region, we have no analytic argument in favour of Fig. [2^, or Fig. [2]d. 

ii) p > 2X n (weak truncation): in this case h > X n . Indeed, again from E[A 2 ] < A„, we deduce p— E[Xlj> 
p — X n > X n > E[Xl]. Here, the correct profile of \x n — E[A, 2 ] | is certainly the one depicted in Fig. ^k. 

Hi) X n < p < 2A„: in this case it is difficult to conclude anything about h, except that by continuity there 
exists a value A„ < p*{X/ n )) < 2A„, possibly depending on A( n ), such that h < X n 4^ p < p*. 
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\xl - EK 





VWl] VP 



b) 



Fig. 2 - Qualitative behavior of the function x n i— > \x\ — E[Jf^]| within the support {—y/p, +y / p)- Note 
that, depending on p and E[X^], the function might have one or two maxima. Specifically, plot a) 
represents the function if p — Epf^] > E[X^], while plot b) applies if p — E[X^] < Epf^]. 



To conclude, the estimate obtained from Holder's inequality is certainly as strict as needed for eq. (1.2) 
to hold true only in case of strong truncation, i.e. for p < X n . In addition, there is a crossover region where 
the same estimate might be sufficiently strict, while it becomes definitely too loose in the region of weak 
truncation. 



3.3 Second argument in favour of eq. (1.2) 



In order to extend the above proof to the region < p < 2A n , we work on the integral representations of 
E[X^] = X\{a nn /a) and E[X^] = X n (a n /a). In terms of these, A„ reads 



A„ = 



A; 



- 2 



(3.14) 



Now, we observe that independently of v, a n -k is bounded from above by (p/X n ) k p a n:p for any p < k. 
Indeed, since x £ B v (p) —J~p < x n < -Vp, we have 



O-n-.k 



A „ . , , p k 



JB v (p) 



dVx HT TT S ( X m,X m ) < 

p „ 1 



x 2p v pk-p 

d x — J^J 5(x m , X m ) — j._p tXn.p ■ 

P m=l A » 



Thus, we immediately obtain 

xt. 



A„ < 



Xn 



a 



< , if p < 2A r , 



(3.15) 



(3.16) 



This conclusion is somewhat conservative, as indeed A„ < p < p*, being implicitly defined by the 
non-linear equation p* = 2{A„ + E[X^ | X e B v (p*)}}. By continuity, the latter is certainly fulfilled by some 
2A n < p* < 4A„. The argument presented here does not apply for p > 4A ra . 



4 Weak truncation expansion 

In order to study the variance of the square components of X in the regime of weak truncation, we need 
to develop an appropriate formalism. To start with, we observe that the constraint X € B v (p) becomes 
increasingly unrestrictive as p — > oo. As a consequence, we have the asymptotic factorization 

V 

(v) / \\ P»max J {A J } TT (1) / \ \ / , -, \ 
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Fig. 3 - Slicing B v (p) orthogonally to the n th coordinate axis X n amounts to representing it as an 
uncountable union of Euclidean balls, living in the (v — l)-dimensional subspace orthogonal to X n , 
with square radius varying from to p. 



The larger is p, the less is the error made in approximating Qti;k 1 ...v.k v by its factorized counterpart. We 



aim at characterizing the corrections to cq. (4.1 1 when p is large yet finite. Actually, we are not interested 
in a complete factorization of the Gaussian integrals: if p ^> X n just for some 1 < n < v, we look at the 
partial factorization occurring along the n th direction. Note that: i) in the regime of weak truncation, every 
rational combination of Gaussian integrals — such as A„ — is led by its factorized counterpart; as we shall 
see, the latter is subject to relevant simplifications in case of ratios of integrals; it) 1-dimensional integrals 
cannot be further simplified, as they amount to lower incomplete gamma functions, 

a«(p;A n ) = -^7(* + ^2^) ' j(s,x) = £ dti*-^* . (4.2) 



4.1 Expansion of Gaussian integrals 

In order to present the idea, we first focus on a. If p X n for some 1 < n < v, we slice the integration 
domain orthogonally to the n th direction, as depicted in Fig. [3] From a geometrical point of view, this 
corresponds to representing B v (p) as an uncountable union of (v — 1) -dimensional Euclidean balls, i.e. 

Bv(j>) = U {y e RV : tf» = x ™ > feo e B -^p - *»)} • ( 4 - 3 ) 

Such technique has been first considered by Ruben [3] with the aim of obtaining an integral recurrence 
relationship on the dimensionality of a. Interpreting the integration domain in terms of eq. (4.3) indeed 
yields 

« W (p;A)= [ + ^dx n 6(x n ,X n ) a^- l \p-xl;\ {n) ). (4.4) 
J-Vp 

Since a(p — x\\ A(„)) is a smooth function of its first argument p — x^, we propose to expand it in Taylor 
series around x\ = + , 

00 / ^ % f~ / 2 \ ^ 

cfr-V(p - xl- X {n) ) = p { -=^-X k n (g) d k p a^\p- X {n) ) 

^^Xin^^f^Xtf^t^p-,^), (4-5) 
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with the functions rjk defined by 

1, k = 0, 

vi v \p-A) = { , k h ( 4 - 6 ) 

[aW(p;A)]-V9^)(p;A), k>l. 



When inserted into eq. ( |4.4[ ), the Taylor series turns into a weak truncation expansion of a, namely 

fe=o ' V " / 



= a«(p; A„) a^Cp; A (n) ) - ^ a W(p; A„) a^fa A (n) ) ^(p; A (n) ) 

+ ^««(p;A„)a( t '- 1 )(p;A (n) )^ 1) (p;A ( „ ) ) + O 0f) . (4.7) 

Although a complete factorization into functions of A„ and A/ n ) is not exactly realized at finite p, we see 
that it occurs at each order of the expansion. We warn that eq. (4.7) has been obtained upon bringing an 
infinite sum under an integral sign. Such exchange of limits is delicate, so it is not a priori evident whether 
the resulting expansion converges or approximates its target just as an asymptotic series. We shall come 
back to this point later on. We also stress that, while power counting is performed by factors of (A^T 'p) 



{XJf 

additional powers and exponentially small terms in p are still hidden within the coefficient functions 2 
and rj k v l \ 

To simplify the notation, in the sequel we drop all function arguments, whenever this does not generate 
confusion. Thus, we shorten eq. (4.7 1 to 

a« - a^a^" 1 ) - ^ a^a^ ^ + -% a&> r,^ + O ( . (4.8) 

P 2 p 2 V P J 

The same technique can be straightforwardly applied to a n:p . For instance, we have for p = 1, 2, . . . 

«s° = «s ^ - x - * {v - i] v ( r i} + \ 4 nt i] + o ( 4 ) , (4.9) 

P L P \P J 



Since the above expansions are all based on eq. (4.5), the coefficient functions r\ k arc the same independently 
of p. By contrast, the multiplicity of the index n of the 1-dimensional integrals contributing to each order 
is shifted forward as p increases. 

Now, when it comes to expanding Gaussian integrals with more than one index, the above procedure is 
carried out in a slightly different way. For instance, in order to expand a nm we need to take into account 
factors of x\ and x 2 m under the integral sign. Accordingly, we slice B v (p) sub sequ ently along the n th and 



directions under the assumption p 3> max{A„, A m }. The analogous of eq. (4.4) reads 



a£l(p;A)= / dx n Y l S(x n ,X n ) / dx m -f k 8{x m , A m ) a ( "~ 2) (p - x 2 n - x 2 m \ A (mn) ) , (4.11) 



-Vp a " J-Vp " 



2 The reader will observe that the coefficient function oS v 1 ' showing up in each term of the expansion is totally useless, as 
amplifies with the one attached to r/ k v 1 ' . ! 
considered, as we shall see in sects. 4.2 and 5. 



it simplifies with the one attached to r/ k v 1 ' . Such redundancy is real, yet it turns useful when ratios of Gaussian integrals are 
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with A( nm ) = {Xk}k^n,m- Again, we expand a(p — a; 2 — x^;A(„ m )) in Taylor series around the point 
x 2 n + x 2 n = + , thus obtaining 

a^- 2 \p-x 2 n -x 2 m ;\ {nm) ) = ^t}^ {x l +x 2 m y d J pa (v-2) {p . X{nm]) 



3=0 J 



E ^# E ( l ) < k *%- k H« {v - 2) (p; V»>) 



^-)( P ;A (nm) )£^£ 

j=0 • y ' fe=0 v 7 v p 



(S)*(sr^<«w- <«.*> 



Here we have also used Newton's binomial formula to express each term of the series in products of powers 
of x 2 and x^. Inserting this expression back into eq. (4.11) yields 

= 4 1} - ^ a£> a« aC- 9 ) rf*> - ^ a« at- a > ^" 2 ) + f . (4.13) 

p p \p 2 J 



Eq. (4.12 ) can be also used to obtain alternative expansions of a n:p if p 3> max{A n , A m } for some m, e.g. 
a W = a« aWa^) - ^ a $ a (D a<»- 2 > - -J a« a« a<"- 2 > ^" 2) + O (^j , (4.14) 

^ = aSl a<- a > - ^ «£> a« ^" 2) - ^ a« a« a^ 2 > ^ + O ( *) , (4.15) 

P P \P J 



4.2 Exercise: relative amplitude of variances and covariances 



We make use of the above expansions to qualitatively compare the correlations X 2 has with itself and the 
other square components of X in the regime of weak truncation. This is rather instructive, because it shows 
how analytic cancellations occur in the proposed formalism. In addition, the exercise inspires the following 
unproved 

Conjecture 4.1. If X ~ Af v (0, A) with A = diag(A), the covariance matrix of the vector {X 2 }^ =1 condi- 
tioned to B v (p) is diagonally dominant, i.e. 

var(X 2 | X G B v (p)) > ]T \cov(X 2 ,X 2 m | X G B v (p))\ , p G M+ . (4.16) 

m^n |~7| 



Eq. (4.16) is supported with no exceptions by extensive numerical tests. As observed in ref. [5], it entails the 
inequality v~ l y] k _ 1 pk < p/( v + 2), where p = {pk}k=i denotes the eigenvalue spectrum of the covariance 
matrix &13. 



In terms of Gaussian integrals the scale invariant observables we focus on are 

r c.i / c..i \ 2" 

r$ = ^var(X 2 |Xe^(p)) = ^ 



0) 
a(«) 




rW = ~cov(X 2 n ,X 2 m \X eB v (p)) 



p p 



0) 

«(«) 



(«) O) 



n =/= m . 



(4.17) 
(4.18) 
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icr J 



10"' 



-S io" 




Fig. 4 - A plot of |r„m| vs. p/\z at {Ai,A2,A3} = {1,2,3}. Solid curves correspond to computed 
values: markers have been overimposed just to make the legend clear. 



For illustrative purposes, we show in Fig. [4] a plot of \Tnm\ vs. /5/A3 at v = 3 corresponding to the choice 
{Ai,A2,A3} = {1,2,3}. Both T nn and T nm vanish as p — > 00, yet the former vanishes as 1/p 2 due to the 
chosen normalization, whereas the latter is exponentially damped. 

In both cases, 



We use eqs. (4.8)-(4.10) to work out the expansion of T nn and eqs. (4.13 )-(4.15 ) for T 
in order to expand a~ L we rely on the Taylor formula (1 — x)~ x — 1 + x + a 
r„„ we have 



a 



(v) 



a 



(v) 



nn 



A„ 
P 



x n:3 
v(l) 



a (1) a {1) 

lx n:2 "n 



Vi 







O(ar). Thus, with regard to 



(4.19) 



a 



(v) 



a 



(i) 



whence we obtain 



a {1) a {1) 

"n:2 "n 



y(l) 



nn nn 3 



a 



*n^3 
(1) 



a (1) a (1) 



a 



n 

(i) 



(v-l) 

m 



o 



o 



We see that the leading term of T nn coincides with its 1-dimensional counterpart. Since lim 
(2k — 1)!!, we have linip^oo (p 2 /A^) Tnn = 2, and thus we find again 



p— foo 1 



-,(«) p»a« Aj 



(4.20) 



(4.21) 



„(!) - 



(4.22) 



apart from exponentially small terms in p. Analogously, we have 



a 



(v) 



a {1) a (1) 

^ n 

«(!) ttW 



A, 



"n:2 _ "m 
rv(l) 



v(l) 



n 

v(l) 



(t>-2) 

»7i 



Q 



(1) 
rn:2 

K(l) 



a (1) /a (1) 



(v-2) 







(4.23) 
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(v) («) (1) (1) 

L n m 



Q 



(1) (1) / (1) 

L n:2 / n 



(i) (i) / (i) 

m:2 n I 



(o-2) 



(v-2) n / A 2 



(4.24) 



and 0(A//?)-terms, so we are left with 



When eqs. ( 4.23 )— ( 4.24 1 are put into eq. (4.181, an exact cancellation occurs separately among the 0(1)- 

(4.25) 



A 



pO) p> max{A Tl: A m } A n A m 
nm ' J \ o 

P P VP. 

Eqs. (4.22) and (4.251 reflect the behavior observed in Fig. [4| 



4.3 Asymptotic vanishing of rjk from combinatorial arguments 

In order to make the weak truncation expansion effective, we need to characterize the coefficient functions 
r]k and provide an algorithmic recipe for their computation. A trivial property, i.e. the vanishing of rjk as 



p — > oo, can be proved from purely combinatorial arguments based on eq. (2.3). This kind of proof nicely 



follows as a simple application of the scaling eq. (2.1), yet it gives no clue to the vanishing rate of We 
begin with the following 

Lemma 4.1. Given a set of m > distinct indices {ni, . . . , n m } and a corresponding set of strictly positive 
multiplicities {fci, . . . , k m }, we have 



lim pd p o%l k k (p;X) =0. 

p— >OQ 11 mm 



Proof. From eq. (2.3), it follows 



lim p dpCi 



P n\\k\ . . . rimikn 



-Y lim a (v \ , + -(v + 2k) lim afc\ 

m:ki...n m :k m n 9^ ' p^oo ni:kl 



2 * — ' p— >oo 
n— 1 



(4.26) 



(4.27) 



with k = ^2j—i kj. The sum index n in eq. (4.27) can either match one of the indices m, ■ ■ ■ ,n m or none of 
them. Moreover, we have 



lim 



lim a, 

p— 7-00 



(v) 



if n 7^ rij V j = 1 , . . . , m ; 



(2kj + 1) lim qA l. H k if n = rtj for some j : 

p— fOO L L mm 



(4.28) 



as a consequence of the exact factorization of the Gaussian integrals as p — > 00 and the standard formula 
E[z 2n ] = (2n - 1)!! (n > 0), valid for z ~ W(0, 1). Therefore, 



v — m 1 
2 2 



m - 

£(2fc i + l) + -( W + 2fc) 



lim 

p^-oo 



ni :fei . . . rimikr, 



3 = 1 

m + 2k m + 2k 



= 0. 



(4.29) 
□ 
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From Lemma 4.1 we easily derive 

Proposition 4.1. As p — ¥ oo all the coefficient functions rj k vanish. 

Proof. Let us define f k = (pd p ) k a^ and x k — J2n ± m.= i a n}...n k - We first prove by induction that f k is a 
homogeneous linear function of xq, . . . ,x k - From eq. (2.3), evaluated at k\ = . . . = k n = 0, we have indeed 



1 v v 1 

pd p a (v) = — a [v) - - ^ a i V) = n x o - = fi (x , x t ) 



fc=i 



Now, suppose that fk—x is a homogeneous linear function of xq, . . . , x k -x- Then, 
fk = {pd p ) k a^ = (pd^pd^a^ = pd p f k ^(x , . . .,x fc _ x ) 



(4.30) 



= fk-i(pd p x , pdpXk-i) = fk-i pdpt 



P d P a nl..n k ., 



n 1 ...n k _ 1 = l 



The inductive step follows from eq. (2.3 1 and the assumed linearity of f k -i- Hence, we have 



lim (pd p ) k a (v) = lim f k (x , ■ ■■,x k ) 



p— >oo 



p— >oo 



(4.31) 



/fe-i \\m pd p a,..., V Km pd p a^ 



m,...,n k -i 



p^-oo 



= A_!(o,...,o) = o. 



(4.32) 



where the last equality is again a consequence of the homogeneous linearity of f k —i and the second-to-last 
one follows from Lemma 4.1. In addition, we known that (see for instance exercise 13, chap. 6 of ref. [3]) 



3=1 



pi ' 



/ith the symbols [ .] denoting unsigned Stirling numbers of the first kind. Hence, we conclude 

k rn k 



lim rj{^ = lim — ^^(-1) 



p— > oo 



k—j 



3 = 1 



3=1 



lim (pd p ) j a (v) = . 

p— >oo 



(4.33) 

(4.34) 
□ 



4.4 Gaussian representation of r\ k 

The above discussion suggests a convenient way to compute the coefficient functions. We have just seen 
that rj k is a linear combination of / 1; . . . , f k . Moreover, V j > fj is a linear combination of xq, . . . ,Xj. 
We conclude that rj k itself can be represented as a linear combination of xq, . . . , x k . Since we know how 
to compute Gaussian integrals with controlled uncertainty, we have a complete recipe for r/ k , provided we 



determine the coefficients of such linear combinations. To this aim, we concentrate first on the f k 's. Eq. (4.30 ) 
gives the analytic expression of f\. By direct calculation we can also derive the expressions 

v 2 v + 1 1 

h(xo > x u x 2 ) = —x — x 1 + -x 2 , (4.35) 

v 3 3v 2 + 6v + 4 3v + 6 1 
fo{x Q , xi,x 2 , x 3 ) = —x xi H — x 2 - ~x 3 , (4.36) 



A generalization is provided by the following 
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Proposition 4.2. For k > 1, we have 



fk(xo, ...,Xk) = ^d M {v)x e , 



(4.37) 



i=0 

where the coefficients dke(v) are defined by 



dki{v) 



E' : .MlTi- /■■ = » ' 



2' 



(4.38) 



0. 



otherwise , 



(4.39) 



v!! / 1 fc = 

(v-2k)U \ ( w _2Jb + 2)^_i fc>l 

and i/ie symbols { ^} denote Stirling numbers of the second kind. 

Proof. The proof is by induction. We first note that dio(v) = v/2 and dn(v) = —1/2. Hence, for k = 1 
eq. (4.37) agrees with eq. (4.30). Now, suppose that fk is well represented by eq. (4.37) with du(v) and 
as in eqs. (4.38) and (4.39). Then, from eq. (2.3) it follows 



f k+ i = ( P d p ) k+1 a^ = d M {v){pd p )xt = £ dki(v) ]T pd p al 



(«) 



.hi 



£=0 



= X>(«) E !^+2%£U-^ E 4:u 

£=0 fe!...fef=l [ fcf+i = l 

= £(?+'W«) E «£U-?E d «(«) E 



(v) 

ki...ke+i 



«=0 
k 



k\...kg=\ 



1=0 



fei...fc^+i=l 



e=o L 



■fc(i-i) 



(4.40) 



The argument is complete provided we are able to show that dki(v) fulfills the recurrence 

d(k+i)e{v) = (| + i) d k i{v) - ~d H t-i){v) . (4.41) 
To this aim, it is sufficient to make use of the basic recursive formulae {"^J 1 } = m { m } + { m ^i} an< ^ 



(n+V\ 



(") 



ml \m- 



■J. We detail the algebra for the sake of completeness: 
fe+i 

d(k+i)i(v) = E 



t=i 
fc+i 



(-1)' (k + Uft 

-*-M t 



= 2^ vt-e t 



k+l 

E 

t=e 



2' 

(-11 
2* 



t \e 



fcl ft 



2^ 2* 

t=i 



t-l\\i 



+ E 

i=£-l 



k\ ft+i 

£ 
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2' 



( — 1) 

t=i 

k 

E 



2 t+i 



H-t+i 



k) ( t 



2* 



-/t 



E 2(«-^)]^-i 



fcl /A 1 



t \\l 2 



4(£-i)(w) = ( - +1) d ke (v) - -4(^-i)(w) • 



(4.42) 
□ 



In view of eq. (4.33), it is no surprise that the coefficients dke(v) embody Stirling numbers of the second 
kind. Recall indeed that Stirling numbers of the first and second kind are related to each others by the 
inversion identity 



<{j\fc} 



E (-u 



t—k 



Ojk- 



(4.43) 



From Proposition 4.2 and eq. (4.43), it follows 
Proposition 4.3. For k > 1, we have 



^W) = E>^) E 



e=o k 
with the coefficients Cke(v) defined as 



(4.44) 



(-1)' 



,!! 



2 k [v -2{k -£)}V. \£ 
0, 



< £ < k, 
otherwise . 



(4.45) 



Proof. We have all the necessary ingredients to carry out the proof. Again, we detail the algebra for the 
reader's convenience: 

1=1 L J 1=1 L J m=0 



fcl L J m=0t=m 



t\\m 



m— t=m 

L v v 

M Z^ Z^ 



('=1 



1 1 V rn 



a (v) ^ ^ 2 * -&-™M m ) a: »» - a („) ^ _>-' 

rn— t—m v 7 m— 



^ — ^ — 0fc-m 



i y' (-i) r 

aW ^ 2 fc 

m— 



(4.46) 
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Note that on the second line above we could extend the upper bound of the sums over m and t from i to oo 
thanks to the property {^} = if b > a. This in turn allowed us to perform the sum exchange on the third 
line. By a similar argument, on the last line we could reduce the upper bound of the sum over m from oo 
to k. □ 



Obviously, computing eq. (4.44) becomes increasingly demanding for larger values of k. Nonetheless, many 



contributions to the sum on the r.h.s. coincide. In particular, all Gaussian integrals with the same index 



multiplicities contribute equally, thus we can recast eq. (4.44 1 to the computationally cheaper expression 



z — ' \mi, . . 

mi . . .m v =0 



,mi-\-...-\-7n v " 



l:mi...v:m v 



(p;A) 



(»)(p;A) 



(4.47) 



4.5 Asymptotic sign of 77^ 

In this and next subsection we put the combinatorial approach on hold and work on the integral representation 
of the coefficient functions. A first property which turns out to be essential to the last part of the paper 
concerns the sign assumed by rjk as p ~ > 00. In regard to this, we state the following 



Proposition 4.4. As p — ¥ 00 the sign of r]k becomes independent of v and A. In particular, we have 



hm sign^ ) (p;A) = (-l) fe - 1 

p— >O0 



(4.48) 



Proof. We first express a in spherical coordinates, i.e. we perform the change of integration variable x = ru 

: \\z\\ — 1} (in the sequel we write 



in eq. JL4| ), with r = \\x\\, u € dB v {l) and dB v (l) = {z e M? 
d v x = r v ~ 1 drdu; here du embodies the angular part of the spherical Jacobian and the differentials of (v — 1) 
angles). Thus, we have 



« w (p;A) 



1 



2I>/2)|A| 1 /2 



dr r^ 1 



exp 



2 V(t 



V(u) = it T A 



u, 



with M representing the uniform average operator on dB v (l), namely 
»r„i _ I>/2) 



27T"/2 



du g(u) . 



(4.49) 



(4.50) 



3B„(1) 



In order to compute rjk, we differentiate a under the integral sign. The first derivative evaluates the ra- 
dial integral at its upper limit, while the remaining k — 1 ones distribute according to the chain rule of 
differentiation. Explicitly, we have 



p k d k p a^{p-\) 



2'"/ 2 I>/2)|A| 1 /2 

(_l)fc-ip«/a 
2«/ 2 I>/2)|A| 1 /2 

(-i) fc -y/ 2 

2"/ 2 r(t»/2)|A| 1 /2 



Qk-l V/2-1 



exp 



k-l 



Q 



k-l 



P V{u) 



exp 



k-i- 



P V{u) 



exp 



pV(u) 



f>=v/2-l 



Here x n = x{x + 1) . . . (x + n — 1) denotes the n raising factorial of x and 

k 



Q k (x,a) = ^2 



1=0 



a e x k - e 



6=ti/2-l 

(4.51) 
(4.52) 
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is a polynomial in x of fc th degree, differing from a Newton polynomial for the presence of a 1 in place of a 1 
We note that the coefficient of the leading term of Qk(x,a) is ( )a>° = 1. It follows that 



p k d k p a^{p-\) 



p-^-oo 



(-1) 



k-l 



r(</» + i)|A|i/2 V2 



;) 14 P{u) k - 1 e X pl [ -'^^} (4.VS; 
Eq. (4.481 follows from the positiveness of V(u). □ 

maxfc{Afe}, we obtain as a by-product an estimate of the exponential 



pV(u) 
2 



Since V(u) > AjJ^, being A n 
damping of r)k , namely 



/d>^(p;A) 



W 2 



< 



2V2r( u /2)|A| 1 /2 



M 






C -P/2A 






(j>=v/2—l 





(4.54) 



4.6 A convergence estimate for the expansion 

We come back to the issue raised at the beginning of this section: is the weak truncation expansion con- 
vergent? It is not difficult to see that the answer lies specifically on the behavior of r/k as a function of k. 
Eq. (4.54) shows that the relevant information is brought by Qk-i(pV(u)/2, —<j>), particularly by its relative 
minima/maxima. As k increases, the position of the latter shifts towards larger and larger values of p, while 
their absolute size increases. In other words, however we choose a reference scale p > we always find k 
such that argmax{|ry fe ( j o; A)|} > p and max{|?7 fc (p; A)|} > max{|?7^(p; A)|} for k > k. For this reason, the 
convergence issue reduces to quantify the increase rate of r)k as a function of k. 

More quantitatively, we first notice the inequality 7(0, x) < x a / a. In order to prove this, we observe that 
a convenient representation of the lower incomplete gamma function is provided by (see for instance sect. 6 
of rcf. [7]) 

7(0, a) = -x a e~ x M(l,l + a,x), (4.55) 
a 

where 

M(a,b,x) = iF^a-Mx) = V (4.56) 

is the confluent hypergeometric function originally introduced by Kummer. Since 1™ = n! and (1 + a) n — 
(1 + o)(2 + a) . . . (n + a) > n\ if a > 0, it follows 

M(l, 1 + a,x) = y - ^— ^- = y - xU <Y^-= e x . (4.57) 
v ; ^ (l + a) n n \ ^ (l + a) n ^ n\ ' 



We can use eq. (4.2) and the above inequality to establish an upper bound to the 1-dimensional Gaussian 
integrals, namely 

1 / \ fc+1/2 



Now, let us denote by X n -.k the weak truncation expansion of a n: k, i-e. 



x i v .l (p; A ) = a(t,-1) (p; A («)) 53 ( v - ) "tlwOnKitf-VfaXw). (4-59) 
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In view of eq. (4.58), an absolute estimate to X n -k is given by 

p 



I^(p;A)I <E 



p=0 



p\\p 



X n)a^( P ; A (n) ) Vp v - 1] (p; A (n) ) 



< 



1 



2tt \\ 



fc+1/2 oo 



' p! (» ■ 



(4.60) 



From the above inequality we see that a less than factorial growth of r\ v with p would make the r.h.s. of 
eq. (4.60) convergent. In order to estimate \p p d^a\, we make use of eq. (4.51). Since dB v (l) is a compact 



domain, we can get rid of the angular average by defining 



u = arg max 

u£9B»(l) 



Qp- 



whence it follows 



P T(u) 



v/2 



< 



2"/ 2 r(?;/2)|A| 1 /2 



exp 



ip-i 



P V(u) 



P V[u*) 



exp 



pV(u*) 



(4.61) 



(4.62) 



We have already noticed that V(u) > X n ^ x for all u 6 dB v (l). If we also consider that IAI 1 / 2 > A^, being 
Ami n = m infcjAfc}, then we have [T^u)"/ 2 ^! 1 ' 2 ] -1 < (Amax/Amin) 1 / 2 . Multiplying and dividing the r.h.s. of 
eq. (|4.62| by V{u*) v ^ 2 leads to 



^^"'(ffA) 



< 



< 



r(«/2) 



Qp-i 



pViu* 



exp 



pV{u*) 



1 



A n 



r(«/2) V K 

Accordingly, we obtain the estimate 



v/2 



max 

x e 



< {^IQp-i^-^le-} 



1 



An 



(«-l)/2 



- 1)/2) \A min ; 

where we have set A m ; n = minj^„{Aj}, A max = maxj^„{Aj} and 
1 



p >^ fc+1 / 2 ~ C<">(p) 



C {v Hp) = - max 



(4.63) 



(4.64) 



(4.65) 



^*=(t<-3)/2 



It will be observed that in order to arrive at eq. (4.64), we have gone through a rather long inequality chain, 



so it is not clear whether the resulting upper bound is finite. For the sum on the r.h.s. to be convergent, it 
suffices that 3 e > 0, A > and po such that C(p) < Ap~ e = m(p; A, e) for p > p . If this holds true, then we 
have C(p)/p < ^C(l + e ) + const., where £ denotes the Riemann zeta function. In order to evaluate 

C{p) analytically, we need to solve a polynomial equation of degree (p — <f>*) for p » 1. An alternative 
approach is to compute C(p) numerically for a set of sufficiently large values of p and then try to fit data 
to a model such as m(p;A, e), with the parameters A and e depending in general on v. In Fig. [5] (right) 
we plot numerical determinations of C(p) for v = 2, . . . , 6. We observe that C(p) is monotonic decreasing 
for v < 5 and monotonic increasing for v > 6. In Fig. [5] (left) we report the values of the fitted parameters 
A and e together with the corresponding x 2_ vahies (the range chosen for all fits is p £ [50,100]). Our 
numerical experiments suggest that the weak truncation expansion converges uniformly in p at least for 
v < 5. Nevertheless, the argument is not conclusive, since it assumes that we can extrapolate the fitting 
model to p — > oo, which is not mathematically rigorous... 
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Fig. 5 - Numerical computation of C' v '(p) for v — 2, ... ,6. We fit each curve on the right to a model 
m(p; A, e) = Ap~ e , with A and e depending in general on v. Fitted parameters are reported on the 
left, together with the corresponding x 2_va lues. 



5 Variance reduction in the regime of weak truncation 



If we look at eq. (4.21), we see that the next -to-leading contribution to the rescaled variance T nn is the 
sum of a few ratios of 1-dimensional Gaussian integrals, all proportional to rji- The subsequent terms of 
the expansion have an increasingly complex structure. Each power of X n /p couples to several products of 
coefficient functions always combined so as to give the correct power counting. If f(p;X) is a generic 
observable, its weak truncation expansion reads 

oo , . v q 

/(p;A)=$>iW^j J2 ^Hp,vvt 1) (p^(n ) r°...vt 1) (p-,\n ) y«- (5.1) 

9=0 ^ P ' eSS? 

We denote by e = {e , . . . , e q } the exponents of 770, . . . , r] q and by S™ the set of all possible e's corresponding 
to an overall power counting m, namely 

S q n = {eeK +1 : V q {e)=m) , (5.2) 

with the power counting function V q (e) defined as 

V q (e) = Y,ke k . (5.3) 

k=l 

Note that if m < q and e G 5™, then e m +i = . . . = e q = 0. In this case, we interpret as 2^ m '^ eo ' " ,em ^. 

Recall also that since r)o = 1, eo never contributes to the power counting. For later convenience, we define 
e = {ei, . . . , e q }. If e € 5™, with abuse of notation we also write e € <S™. Strictly speaking, the presence of 
170 in eq. ( |5.1| ) is necessary in order to properly take into account the leading order of the expansion, namely 

00 

hm f(p;X)= Km ^4 0;{eo}) (p;A). (5.4) 

e =0 



The sum over eo in eq. (5.1) extends in principle from to 00. We use eo to enumerate all contributions 
proportional to rf^ . . . rf q . Accordingly, the information concerning the maximum value taken by eo is hidden 
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within 2y?' s ' 
cases of V nn 



We finally stress that Sj 3 '-- 1 depends in general upon all the components of A, yet in the specific 
and A„ (which are the ones we are interested in) it depends only upon A„. 



Suppose now that /(p; A) and g(p; A) are two observables, which we expand according to eq. (5.1). It is 
not difficult 
/ + g and / 



to prove that the convolution rules needed to obtain the expansion of the algebraic combinations 
g are similar to the Fourier's ones. Specifically, we have 



(5.5) 



and 



;c) 77 ( m; d) 

a 



(5.6) 



\rn=0 cG5« dGS™ 



where 6 CJ i = Y[1 = o3c k .d k is a vector generalization of the Kronecker symbol. Eqs. (5.5| and (5.6) are 



sufficient to derive the weak truncation expansion of A„ to all orders and to consequently prove our main 
result, represented by the following 



Theorem 5.1. As p — > oo, the sign of the coefficient function E^'~\p; A n ) is given by 
1 nfes) 



lim^f(p;A„) = (-l)£U 

) — text ' " n 



p— too 



(5.7) 



This, in conjunction with eq. (J i .J i 8), implies that all terms of the weak truncation expansion 



of A n become negative at sufficiently large p. 



Proof. We proceed in subsequent steps. First of all, we observe that A n can be written as 

J 3 n = a$ a<»> - ^ai"' - 2 ^aW , 

To work out J) n , we first review the expansion of a n -k- Explicitly, we have 



A„(p;A) = ^Dn-D 



9=0 



(fc+9) a ^9 



(5.8) 



E 



9=0 

Hence, we deduce 



n:(k+q) 



e (n^,oUe„ 1 [r ) r...[<- i) ' 



1 



,.(!) 



This allows us to derive the expansion of the product ot n:r a n - a . From eq. (5.6) it follows 

9 



1— '«« !r Q!n: 8 



E 



<W+m (1) (1) , 

l\m). n '-( r + i ) n:(s+m)t 



„(«-l)l2 



e e ^ (n^°) ( ff^.o) 

cS5* dG5™ \i=0 I \ fc=0 / 



(5.9) 



(5.10) 



(5.11) 
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The inner sums can be performed exactly. Indeed, thanks to the Kroncckcr symbols 5 Ct i and Sd m ,i, non- 
vanishing contributions group according to whether i = m or £ ^ m, namely 



5, 



„(i) 



£,m=0 



l\m\ n: ( r + e ) n:{s+m) 



0...q 0...q 



(5.12) 



We see that the Kronecker symbol 5 q ^+ m makes both groups of terms vanish unless V q (e) = q, as intuitively 
understood. Conversely the only elements e S which result in a non-vanishing coefficient function 

si^~' Q „. 3 are either those where two different exponents eg, e m equal one (with £ + m — q) while the others 
vanish, or those where e q /2 — 2 and e, = for i ^ q/2 (of course the latter contribute only when q is even). 
From eq. (5.12|, we immediately obtain 



£,m=0 



> a (1) -a (1) a (1) 

(<+2) "n:m u ra:(£+l) u n:(m+l) 



a (1) 



a...q 0...q 



(5.13) 



The above expression depends upon p essentially via the integrals in parentheses (the overall factor [cr - " -1 )] 2 
is irrelevant to our aims). Since a n - r — > (2r — 1)!! as p — > oo, we have (a„ : (£ + 2) a n - m — a n -a+i)a n: ( m +±) — 
2ck„ : (£+i)Q!„: m ) — > (£ — rn)(2£ + 1)(2£ — l)!!(2m — 1)!!. In particular, this quantity vanishes for £ = m, thus 
making the first term in square brackets never contribute as p — > oo. A little additional algebra yields 



hm S^ } (p;A) =4V V <5 9/+m (^-m) 2 ^ IT <S ei , . 

^=0 m^O i^£,m 



(5.14) 



We notice that the r./i.s. of eq. (5.14) vanishes always for eo > 2, but not necessarily for eo = or eo = 1. 

As a second step, we work out S)d- To this aim, we first need to evaluate 3 As already done in 
sect. 4, we make use of the Taylor series (1 + x)^ 1 = J2^La(~^) P xP ■ F rom 

(-i)g A B y4 M ; 



it follows 



9=1 



p J a 



(i) 



(5.15) 



fa^o^- 1 )!- 1 



E(-d p E 

p=0 £ 1 ...£ p = l 



(-1) 1 p ( K\ "n:li ■ ' ' a n:£ p (v-1) 



(v-1) 



A, 



q q 



E(-!) p E w.. + ^n 

P =o e-L...e p =i 3=1 



1 "n-Jj (v-1) 

£A aW % 



(5.16) 



On the second line we have reduced the upper limit of the sum over p from oo to q. The reason is that 
8q,ei+. ..+tp — f° r 1 < Pi owing to l\ + . . . + £ p > p. On expanding the sums over £\, . . . , £ p , we see that all 
terms proportional to rf^ . . .r\q q for some e coincide. Since permutations of £x,...,£ p corresponding to the 
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same e give all the same contribution, the latter turns out to be multiplied by an overall numerical factor 
which is at most p\. Of course, permutations of equal indices contribute only once. Therefore, a correct 
counting of that factor amounts to the multinomial coefficient p\/{e\ \ . . . e q \), with the constraint ^2 e j = P- 
In other words, we have 



9=0 v p ' eeS « \e ,...,e fc / j=Q 



1 a. 



(i) 



and consequently 

^ = (-1)^^,0 



eot • • • i efc 



n 

3=0 L 



1 a 



(1) 

n:j 



Now, we obtain by convolving eq. (5.18) with itself. This yields 



(-i) a =' ei <s eo ,o n 



1 a 



(1) 

n:j 



with the coefficient ^(p ;e ) defined by 



From eq. (5.19), it follows 



hm H ^ = (-l)EL 1 e^ "Q 



p— >oo 



(2j - 1)!! 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



As a third step, we convolve eqs. (5.14) and (5.21). In this way we obtain E^'~' directly in the limit 



p — » oo. The algebra is just a little bit intricate, so we detail it. First of all, 



lim ~>lf = 4 ^ <W™ E E ^q.co* 



£,m=0 



f r-1 



o..e 



(_i)Er =1 *. JJ 
fe=i 



(2fc- l)!P dfc 



,(2r-l)!! (2s- 1)!! 



A'! 



^(m;d) 



(5.22) 



Owing to the Kronecker symbols, we pay no price if we introduce a factor of (— l)5D*>=o Cfc = 1 within square 
brackets. For the same reason, we can also insert additional factors of [(2k — l)!!/fc!] Cfc = 1 for k ^ r,s 
without pay. Hence, 

q r (2k- l)!P £fc q 
k\ 



hm = 4(-l)£*=i e * TT 

i — ±r>D n J- 



p— fOO 



fe=l 



E w. E E (-^"•We 



e,m=o 



o...e 



m\d) 



=4(-i)^^ efc ^ ( ,), 9 n 



i^r,s 

(2k- 1)!! 



n c fc g 



r-l 



fc! 

0..J 



E <w E (-i) 00 ^ 



EE^.''+ s ( r_S ) 2(5c - l(5c - 1 n ^ 4> 0* (r ' 



(5.23) 



r=0 s=0 
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with #e.c = Yii=i ®ei,£i being a vector generalization of the Heaviside function 



,b — 



1 if a > b , 
otherwise . 



(5.24) 



From eq. (5.23) it follows 



p— foo 



e =0 



k=l 



(2k - 1)!! 

jfei 



(5.25) 



with 



-1 r-1 



£,m— r— 1 s— 1 



,e*,...,e a }) 



and 



(g;e) 



{e 1 ,...,e^-l,...,e a }) 



(5.26) 



(5.27) 



«,m=0 



(5.28) 



In order to complete the proof, we need to show that Q^ 1 '^ < fl'f'^ Ve £ 5|. To this aim, we hnd it 
convenient to use a slightly different representation of <J/(P ;e ), viz. 

* (p;5) =<wee ( ELiCk )( ^ {ek - ck) y erc . 

n=lciS" V e i -ci,...e, -c g y 

Suppose that q > 1 for some 1 < ^ < q. Then, 

1 Vj/('? _£ ^ e i'---' e '!- 1 '---' e 'i}) 

= ^(s]jEEf^ =lC(: )f - c fc) - 1 

yCl , ■ ■ • , Cq/ \ei Cl , . . . , e n C n 1, . . . Bq Cq 

On performing the change of variable ci cg_ + 1, we easily arrive at 

< W *(«-*;{e 1) ...,e,- 1) ...,eJ) = g £ £ «_ <* ^ ( ELl(^ " Cfc) 

f^cts* ELi c * \ci,---,c 9 y \ei -ci,...,e,-c, 
Analogously, for e s > 1, e r > 1 and l<s<r<g, we have 

6> es i6» e 1 ^(9-s-n{ei,...,e a -l 1 ..., er -l,...,eJ) 



l...g 

?e„ ,c„ + l J J ,< 



(5.29) 



(5.30) 



3 -P,(s), g 



E E 



(ELi c *;) (ELi c fe - 1) V c l> • • • > <V \ e i - Ci, . . . , e 9 - cj e ' c E '=! Cfc ' 2 



(5.31) 



Note that the lower limit of the sum over t in cq. (5.30) can be reduced from £ to 1, since eg = if e G <S* 
and t < I. The same cannot be done in eq. ( |5.31 ) without increasing the resulting sum. Therefore, 



(g;e)_, \ - ELi ^ f ELi c k\f El=i( e k-c k ) 

1 " ^' q hU TU^\c 1 ,...,c q ){e 1 -c 1 ,...,e q -c 



(5.32) 
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and 

-1 r^r-l 



n (q;e) ^ r ST^ 2-^1=1 l^r=l l^s=l °i,r+s{r - S) c s c r 

9CeU hh (EJ=i«*)(ELi^-D 



c fc \/ T,t=i( e k - c k ) 
ci, • ■ ■ , c q J \ei - ci, . . . , e 9 - 

Now, it is immediate to prove that 



(5.33) 



q l-l i — 1 g / g \ 

EEE^( r - s ) 2c - c '' <E^ C M E^- 1 ' Vce5^. (5.34) 

1=1 r=l s=l 1=1 \k=l / 

Indeed, given c G <S* each non-vanishing contribution c s c r > with s < r is weighted by (r — s) 2 on 
the l.h.s. and by (r 2 + s 2 ) on the r./i.s. The remaining terms on the r.h.s. are Y^e=i ^ 2c t{ c t ~ 1) — an( i 
( ? 2c <3(Sfc=i c fc ~ 1) — 0- This concludes the proof. □ 



6 Concluding remarks 

Conditioning a vector X ~ M v (0, A) with A = diag(A) to a centered Euclidean ball B v (p) of square radius p 
affects non-trivially the covariance matrix of its square components. Since the conditional moments of X can- 
not be calculated in closed-form, the only viable approach (besides numerical computation) to characterizing 
the truncational effects consists in establishing analytic bounds to the conditional correlations (variances and 
covariances) of the square components of X. Such estimates are also referred to in the literature as square 
correlation inequalities. 

In this paper, we specifically focused on the conditional variances. In particular, our aim was proving 



eq. (1.2). The analyses presented in the previous sections go in this direction, yet they do not solve the 



problem in a conclusive way. The arguments proposed apply in the opposite regimes of strong and weak 



truncation. For < p < 2A„, eq. (1.2) is easily proved. A bigger effort is required for p ^> X n . Nothing 
is said regarding the intermediate region. We conclude with two major criticisms, representing at the same 
time an outlook of future research: 

• the weak truncation region is not sharply defined: the asymptotic property stated by Theorem 5.1 is 
certainly sufficient to prove that the p th order of the weak truncation expansion of A„ is negative at 
p > p* for some p*, but the theorem does not provide any estimate of p*. A better characterization 
of the coefficient functions % and 5a„ far from the asymptotic regime would help identify precise 



conditions to extend the proof of eq. (1.2) to large yet finite values of p along the same lines of 
Theorem 5.1; 

we also lack a general proof of convergence of the weak truncation expansion. The argument presented 
in sect. 4 suggests uniform convergence in v < 5 dimensions, but it is based on a numerical estimate 
of the vanishing rate of the p th term of the expansion, which cannot be legitimately extrapolated to 
p — > oo. 

The weak truncation expansion of a given observable / (built from Gaussian integrals otktm...) is to all extents 
a perturbative expansion around the factorized value / takes as p — > oo. As such, it is affected by the usual 
problems encountered with perturbative expansions. Having proved a property of A„ to all orders represents 
the main (non-trivial) contribution of the present paper. 
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